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We calculate the entropy of a two-dimensional Fermi Liquid(FL) using a model 
with a contact interaction between fermions. We find that there are T 2 contributions 
to the entropy from interactions separate from those due to the collective modes. 
These T 2 contributions arise from non-analytic corrections to the real part of the self- 
energy which may be calculated from the leading log dependence of the imaginary 
part of the self-energy through the Kramers-Kronig relation. We find no evidence of 
a breakdown in Fermi Liquid theory in 2D and conclude that FL in 2D are similar 
to 3D FL's. 
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The unusual nature of the normal state properties of the high temperature supercon- 
ductors (HTS) has generated a new interest in the metallic phase of strongly correlated 
electronic materials. In particular, much attention has been focused on the existence |1|-|3| 
or non-existence @-§| of a Fermi liquid phase for these systems in two dimensions (2D). 
This controversy as to the existence of 2D Fermi liquid(FL) is motivated by the difficulty 
of fitting some experimental data on the HTS materials with conventional FL expressions 
and also by the property of one dimensional systems that the ground state of a system of 
interacting fermions is a Luttinger Liquid (LL) rather than a FL. In particular this has 
lead to the development of the marginal Fermi Liquid (MFL) phenomenology ||[7| which 
has been used extensively to fit data. However there is no microscopic calculation as yet 
which leads to a an MFL ground state. 

The stability of a FL ground state has been studied extensively in the dilute limit and 
for weak coupling || in 2D. In the dilute limit it is possible to show that the particle-particle 
channel diagrams contribute in leading order. The presence of a two hole bound state in 
this channel lead to speculations that this was a possible source of the breakdown of the 
Fermi liquid phase. |TD[ It appears, though, that this bound state only gives rise to higher 
order corrections to the properties of the FL. In weak coupling away from half filling this 
stability of the FL phase of the 2D Hubbard model was also observed in the propagator 
renormalized fluctuation exchange approximation of Serene and Hess. @ In this approach 
all of the known instabilities, superconductivity, spin and charge density waves, and the two 
hole bound state could occur. No evidence for a breakdown of the FL phase was observed 
at quarter filling. 

From these studies we see that the FL phase of the Hubbard model is stable against 
particle-hole or particle-particle fluctuations away from half-filling. From the two-hole bound 
state it was shown [[!]] that this contributed a term of the order |e| 5 ^ 2 to the imaginary part 
of the self-energy, E (p,e), and from Kramers- Kronig a similar term is found in Re S(p,e). 
The two hole bound states are predominantly short wavelength fluctuations, in this Letter 
we investigate the long wavelength fluctuations. We find that they give rise to lower order 
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non-analytic corrections to FL theory than does the two hole bound state. In particular, 
we find that Re 5£(p,£ p ) oc sgn (£ p )£ p 2 and that this term gives rise to a T 2 correction to 
the specific heat, Cy. What we learn from our work and that of references 1 and 2 is that 
the breakdown of FL theory must be more subtle than is found in any of the traditional 
perturbation theory approaches. The present calculation can not address the issue raised 
by Anderson PJlOj as to the validity of perturbation theory in 2D except to say that there 
is no indication of this from perturbation theory itself. 

Apart from the question of stability of the Fermi liquid our results for the corrections to 
Fermi liquid theory are most surprising. In 3D the specific heat is known to have a T 3 InT 
correction, i.e. Cy ~ 7 T + r 3 £> T 3 InT. In the ID system the breakdown of the FL can be 
seen already in the 2 nd order perturbation theory where the specific heat correction is given 
by 5Cy ~ F\d T InT. (This is a clear signal that perturbation theory does not work since 
this is more important than the linear term). In the 2D case we find that Cy ~ 7 T + T2D 
T 2 +... . One might have expected this to have a T 2 lnT correction by studying the 3D and 
ID behavior, this in fact is not the case. 

In order to determine the leading corrections to a 2D FL due to longwavelength in- 
teractions we consider a system of fermions which interact via a two-body potential as in 
Eq.(l). 

\v\ 2 
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Expanding in the particle-hole channel the effect of the interaction may be considered as 
coming from two independent channels, the symmetric(s) and the antisymmetric(a) channels 
corresponding to no spin exchanged and to spin 1 exchanged. Using the paramagnon model 
first introduced by Doniach and Engelsberg ||11]| , for the q dependence of the interaction, 
the value of the interaction in the symmetric channel, V s , is -1/2 and in the antisymmetric 



channel, V a , is 27, where I > 0. [T2|] I is the the strength of the interaction and multiplied by 
the density of states of one spin is the paramagnon parameter, I. The interaction is cutoff 
at \q\ = q c . In order to compare the properties of FLs in 2D and 3D we first calculate the 



single-particle self-energy to second order in perturbation theory at zero temperature. For 
a 3D FL Blaizot and Friman IlIBl found 
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7mE(p, e) = — - — ^sgn(e){[q c v F - |£ p |]e 2 + -|e| 3 + ....), (2) 



where £ p = (p — pp)vp, Pf is the Fermi momentum, and Vp = Pf/ 771 is the Fermi velocity. 
The real part of E(p, e) is determined from the Kramers-Kronig relation, 



B^ t) = L P r ( 3 ) 
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and is given by 

i2eE(£ f p ) = A 3D C P + B 3D C P 3 In |£ p | + .. (4) 
The £ p 3 /^|£ p | term comes from the |e| 3 term in ImS(p, e) and so it is determined by the 



long- wavelength scattering. In fact it has been shown by Moriya that no £ p 3 ln |£ p | terms 
come from finite q scattering. For a 2D FL we find 

4/2 fe 2 + 2£ P (e - Q)ln{max%, \e\])Q{q c v F - \e\) (5) 
+(e - Q 2 ln(\e - ^\)Q(q c v F - |e - £ p |) + •• • 



I™me) = s9<e)^ mv2 



F 



for £ p > 0. This reduces to the well-known results of Hodges et al. and Bloom [16j for 
ImE(p, e = £ p ). Using the Kramers-Kronig relation again, ReEQo, £ p ) is 

ReE(p, Q = A 2D £ P + B 2D sgn(Q£ 2 + (6) 

Only the first term in ImEfjo, e) contributes to ReE(p, £ p ). The linear terms in ReE(p, £ p ) 
for 2D and 3D are effective mass enhancements, which come from all <f s and depend on the 
q dependence of the interaction. Here we are concerned with the corrections to the effective 
mass enhancement terms. As in the case of B 3 f>, B 2 d is determined by long wavelength 
scattering. Comparing the corrections to the effective mass terms in ReE(j9, £ p ) in 2D and 3D, 
one sees that in 2D the correction is non-analytic and comes from the leading e dependence 
of ImE(p, e) whereas in 3D the correction is analytic but comes from non-analytic terms in 



ImS(p, e). This difference between the leading corrections to the linear £ p dependence in 
2D and 3D is due solely to the different phase space. This may be seen by calculating the 
contribution to the spectrum in 2D using the equation, 

Ae p = £ [l-2f f+ i(p.q) 2 W, (7) 

\q\<Qc 

where is the coefficient of the (p.q) 2 term in the effective quasiparticle interaction. In 3D 
this gives the £,pln\£ p \ dependence. [[L7|] In contrast to ID, where the analogous calculation 
already shows that FL theory has broken down, there is no indication of a breakdown of FL 
theory in 2D to this order in perturbation theory. We now consider the thermodynamics of 
a 2D FL and compare the results with 3D. 

Using the RPA approximation the change in the thermodynamical potential due to in- 
teractions in Eq.(l) is given by 



An = k B T 
where 



3 1 

-(ln[l - Ix(q,uj n )} + Ix(q,uJ n )) + -(ln[l + Ix(q, w n )] - Ix{q,u n )) 
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x(g» = 2^ fe* U y (9) 

fp is the Fermi Dirac distribution function, and u n = 27r(n+l)T are Matsubara frequencies. 
When analytically continued to the real uj axis Afl can be easily broken up into a quasipar- 
ticle contribution, AQ qp , and a contribution from collective modes, AVL co n mo( ies- First we 
consider AQ qp which is given by 

roo di^j 

A*V= E / —n B (uj)(F(q,uj) + Ix"(q,uj)) (10) 

\q]<qc 

F( ^^2 tan [ i-M^) ] + 2 tan [ TTW^) ] ' (U) 
x(M=x!{M + *x!'(q,u), (12) 

and ub{oj) is the Bose distribution function. From this the change in the entropy is 
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The two terms on the right of Eq.(13) involve the temperature dependence of x{Qi u ) which 
is weak when /x is kept constant. 
Calculating AS qp for 2D one finds 



AS qp = 72D /T + T 2D T 2 + 0(T 3 
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where A=s or a, ^ s =l, f a =3, and 
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and A a are the scattering amplitudes in the symmetric (density) and antisymmetric(spin) 
channels, n is the density of particles, and Tp = vppp/2. The T 2 term in AS qp comes 



from the non- analytic term in ReS(p, £ p ) as may be seen from the following argument. [18| 
Consider the entropy of a FL whose spectrum is given by 



£p £p ~\~ Asp 
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where Ae p arises from interactions. Substituting this spectrum into the expression for the 
entropy of a non-interacting Fermi liquid and expanding to linear order in Ae p one finds 

C P Ae p df(e) 



as = y, 



T T de 



(19) 



Assuming that Ae p can be expanded in a power series in one finds 



AS = N(0)k B [°° ^ ] 2 g £ o»r (20) 
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One sees that only odd n terms contribute to AS" and that they lead to series of odd powers 
of T. The presence of the T 2 in AS qp clearly arises from the non-analytic nature of the 
correction to the spectrum in Eq.(6). The terms of 0(T 3 ) and higher are a sum of odd 
powers of temperature. 

Carrying out the calculations in 3D one finds that 

&S qp = l3D /T + T 3D T 3 lnT + 0(T 3 ) (21) 

Another difference between 2D and 3D is that I^d depends on the scattering amplitudes to 
all orders whereas Y 3 d involves only the second and third powers of the scattering amplitudes. 
The collective mode contribution in 2D is 

AS colLmode = T"T 2 + 0(T 4 ) (22) 

where 

and c is the velocity of the collective mode given by 

l + I 

c = v F 24 

VT+21 

The collective mode spectrum does not contain any log dependence on |g| and so there is 
no InT contribution to AS co a, mo d e - In 3D AS co ii mo d e ~ T 3 and does not contribute to the 
T 3 lnT corrections except to change the cutoff of the logarithmic temperature dependence. 

Collecting AS qp and AS co u, mo d e together one sees that AS is a power series in temperature 
in 2D. In particular there are no InT terms in AS which implies that, at least to leading 
order, quasiparticle damping effects do not contribute to thermodynamic properties in this 
approximation where the propagators are unrenormalised. The effect of finite quasiparticle 
lifetimes on the entropy may be estimated with Eq.(23) of ref. |]19[ and is given by 



AS damp = £ /~ de^G(X(p, e)), (25) 



where 



G(X(p,e)) = -tan-\\{p,e)), (26) 

1 + \(p, e) 

and 

JmS(p,e) 

A(p ' e) = "^E(M- (27) 

The functional form of the integrand in Eq. (25) is very complicated and analytic evaluation 
is intractable. So we content ourselves with an estimate. G(X(p, e) is a smooth function of 
X(p, e) which goes as X(p, e) 3 for small X(p, e) and is a constant for large X(p, e). In order to 
get an estimate of ASdamp we assume that G ~ X{p, e) 3 for all values of e. Since there is no 
contribution to the integral for large values of \X(p, e)|, this is clearly an overestimate. With 
these approximations one finds ASdamp ~ T 7 ln 3 T which is higher order in temperature 
than the T 2 corrections found above. Lifetime effects lead to higher order effects in the 
thermodynamics than T 2 so that they are much less important in thermodynamics than a 
calculation of S(p, e) would suggest. 

Since we have used a particle-hole expansion in the symmetric and antisymmetric chan- 
nels our results may be easily extended to Landau's Fermi liquid theory by considering 
Eq.(l) to describe quasiparticles with an effective mass interacting via an effective interac- 
tion fipipi) which can be decomposed into two channels, f s (p, p') and f a {p,i?)- This effective 
interaction is a long wavelength limit of the particle-hole irreducible four-point vertex and so 
describes long wavelength properties. This leads to somewhat more complicated expressions 
when V s and V a are substituted for by f s {p,f?) and f a (p,P^)- The Landau functions are 
functions of the variable s = and may be expressed as coefficients in a series of Legendre 
polynomials in which s is the argument. These coefficients are the Landau parameters. 

The present calculation indicates that a 2D FL is very similar to the 3D case and that any 
breakdown of the FL in 2D has to arise from effects which are more subtle than those which 
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give the leading corrections to FL theory in 3D. The logarithmic dependence in Im£(p, e) 
allow us to keep track of the contribution to the thermodynamic properties from lifetime 
effects. We find that in spite of the £ p 2 /n£ p dependence of the relaxation time in 2D their 
contributions are higher order in T than the contributions to the quasiparticle spectrum. 
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